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ABSTRACT 

We use publicly available galaxy merger trees, obtained applying semi-analytic tech- 
niques to a large high resolution cosmological simulation, to study the environmental 
history of group and cluster galaxies. Our results highlight the existence of an intrin- 
sic history bias which makes the nature versus nurture (as well as the mass versus 
environment) debate inherently ill posed. In particular we show that: (i) surviving 
massive satellites were accreted later than their less massive counterparts, from more 
massive haloes; (ii) the mixing of galaxy populations is incomplete during halo assem- 
bly, which creates a correlation between the time a galaxy becomes satellite and its 
present distance from the parent halo centre. The weakest trends are found for the 
most massive satellites, as a result of efficient dynamical friction and late formation 
times of massive haloes. A large fraction of the most massive group/cluster members 
are accreted onto the main progenitor of the final halo as central galaxies, while about 
half of the galaxies with low and intermediate stellar mass are accreted as satellites. 
Large fractions of group and cluster galaxies (in particular those of low stellar mass) 
have therefore been 'pre-processed' as satellites of groups with mass ~ 10^'^ Mq. To 
quantify the relevance of hierarchical structure growth on the observed environmental 
trends, we have considered observational estimates of the passive galaxy fractions, and 
their variation as a function of halo mass and cluster-centric distance. Comparisons 
with our theoretical predictions require relatively long times (^5 — 7 Gyr) for the 
suppression of star formation in group and cluster satellites. It is unclear how such 
a gentle mode of strangulation can be achieved by simply relaxing the assumption 
of instantaneous stripping of the hot gas reservoir associated with accreting galaxies, 
or if the difficulties encountered by recent galaxy formation models in reproducing 
the observed trends signal a more fundamental problem with the treatment of star 
formation and feedback in these galaxies. 
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1 INTRODUCTION 

The observed properties of galaxies have long been known to 
depend on the 'environment' in which they are located. This 
correlation was quantified by early work showing that galaxy 
clusters are distinguished by lower fractions of star form- 
ing, disc dominated galaxies than regions of 'average' den- 
sity (e.g. iHubble fc HumasonI [l93ll : lOemled [l97i : iDresslerl 
Il980l ). In recent years, observational studies trying to assess 
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the role of environment on galaxy evolution have received 
much impetus from the completion of large spectroscopic 
and photometric surveys at different cosm i c epochs (e.g. 
Kau ffmann et al.l |2004 iBalogh et al. I l2004l : ICucciati et al.l 
i2ood r ICooper et al.l 1200a just to mention a few). Despite 
much effort, however, disentangling the environment(s) and 
related physical processes that are responsible for the ob- 
served trends has proved difficult, and their physical origin 
is still subject of an active debate. 

Much of the argument centres on whether these trends 
are the end product of physical processes coming into play 
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only after galaxies have become part of a 'group' or of a 
'cluster' (the nurture hypothesis) , or whether they are estab- 
lished before these events take place, due to galaxy forma- 
tion proceeding differently in overdense regions (the nature 
hypothesis). In the current standard paradigm for structure 
formation, dark matter collapses into haloes in a bottom- 
up fashion: small objects form first and subsequently merge 
into progressively larger systems. As structure grows, galax- 
ies join more and more massive systems, therefore experienc- 
ing a variety of environments during their lifetimes. In this 
context, the nature-nurture debate appears to be ill posed, as 
these two elements of galaxy evolution are inevitably and 
heavily intertwined. 

As for the nurture scenario, a variety of physical pro- 
cesses might be effective in suppressing star formation and 
affecting the morphology of cluster galaxies. Broadly speak- 
ing, these can be grouped in two big families: (i) inter- 
actions with other cluster members and/or with the clus- 
ter potential, and (ii) interactions with the hot gas that 
permeates massive galaxy systems. The influence of these 
physical processes and their characteristic time-scales have 
been studied u sing detai l ed nu merical simulations (for a re- 
view, see e.g. iDe Lucial I2OIII ). This work has shown, for 
example, that the pressure experienced by galaxies travel- 
ling through a dense intra-cluster medium (ram-pressure) 
can be effective in removing the galaxy interstellar-medium, 
thereby suppressing subsequent episodes of star formation 
IIGunn fc Gottlll97d ). Ram-pressure is expected to be more 
important at the centre of massive clusters, because of the 
large relative velocities and higher densities of the intra- 
cluster medium. Galaxies orbiting in massive clusters also 
experience repeated fast encounters with other cluster mem- 
bers. The cumulative effect of these encounters (harassment) 
can drive a strong internal dynamical response, leading for 
example to the transformation of spiral galaxies in to dwarf 
spheroidals l|Farouki fc ShapirolHoSlI : iMoore et aLlirgQa ). 

As discussed above, however, these physical processes 
should be coupled with a history bias that is integral part of 
the hierarchical structure formation. Until about one decade 
ago, this effect was believed to play a minor role: early nu- 
merical work found no dependence of the clustering of dark 
matter haloes o n their properties, such as concentration or 
formation tim43 /Lemson" fc Kauffmannlll999l : |Percival et al.l 
[2003). Taking advantage of high resolution simulations of 
structure formation, however, recent studies have demon- 
strated that halo properties like concentration, spin, shape, 
and internal angular momentum exhibit clear environmen- 
tal dependencies (e.g. lAvila-Reese et ahl l2005h . Haloes in 
over-dense regions form statistically earlier and merge more 
rapidly than haloes in regions of the Universe of average 
density (|Gao et al.ll2005l : iMaulbetsch et al.ll2007l '). This dif- 
ferential evolution is bound to leave an 'imprint' on the ob- 
servable properties of galaxies that inhabit different regions 
at any cosmic epoch. In this context, a crucial missing in- 
gredient for a correct interpretation of the observed environ- 
mental trends is represented by a detailed characterization 
and quantification of the environmental history of group and 
cluster galaxies. 



^ This is usually defined as the time when half of the final mass 
of the halo is first assembled in a single object. 



A few recent studies have touched this issue. 
iBruggen fc De Lucial l|2008l ) combined semi-analytic mod- 
els of galaxy formation with analytic models for the gas 
distribution in clusters to study the ram-pressure histories 
of present day cluster galaxies. They showed that virtually 
all cluster galaxies suffered episodes of ram-pressure dur- 
ing their life-time, and argued that this physical process 
might have a significant role in shaping the observed prop- 
erties of the entire c luster galaxy population. More recently, 
iBerrier et all (|2009l ') and iMcGee et all (|2009l ') have studied 
the accretion history of galaxies onto clusters with the aim 
to quantify the relevance of pre-processing in galaxy groups. 
These two studies make use of different methods, and reach 
different conclusions. In the following, we will discuss find- 
ings from these studies in more detail, and will compare 
them with our results. 

In th is study, we adopt a n approach similar to that em- 
ployed bv lMcGee et al.l (|2009l ). In particular, we take advan- 
tage of publicly available galaxy merger trees, obtained by 
applying semi-analytic techniques to a large high resolution 
cosmological simulation. These merger trees are analysed in 
order to study the history of the environments that galax- 
ies have experienced during their lifetime, as a function of 
the parent halo mass at present day, and as a function of 
present day galaxy stellar mass. The layout of the paper 
is as follows: in Section [2] we provide a brief description 
of the models used in our study, and describe the method 
and definitions adopted. In Section |3] we characterize the 
history of galaxies in terms of their parent halo mass, and 
define times that should play a relevant role in their evolu- 
tion. In Section|4l we discuss the expected trends due to the 
environmental history for galaxies of different stellar mass. 
We also compare these expectations with observational esti- 
mates of the fraction of red and passive galaxies as a function 
of parent halo mass, galaxy stellar mass, and cluster-centric 
distance. Finally, we discuss our results in Section [S] and 
summarize our conclusions in Section |6] 



2 THE GALAXY FORMATION MODELS 

In this study, we take advantage of the publicly available 
catalogues from the galax y formation model presented in 
iDe Lucia fc Blaizod ||2007 ), and applied to the Millennium 
Simulation ( Springel et al.l |2005| ) . This simulation follows 
N = 2160^ particles of mass 8.6 x 10* H'^Mq within a co- 
moving box of 500 /i~^Mpc on a side. It is based on a ACDM 
model with parameters fim = 0.25, Qh ~ 0.045, h = 0.73, 
f^A = 0.75, n = 1, and erg = 0.9, where the Hubble constant 
is parameterised as //o = 100 /ikms~^Mpc~^. The simu- 
lation outputs were used to construct merger trees of all 
gravitationally self-bound dark matter subhaloes down to 20 

0- 



particles, which corresponds to a mass of 1.7x10^° /i-^M, 
These merger trees represent the bas ic input needed for 
the se mi-analytic model described in iDe Lucia fc Blaizod 
|2003). 

Down to the resolution limit of the Millennium Simu- 
lation, this model has been shown to provide a reasonable 
agreement with a large variety of observational data, both 
in the local Universe and at high redshift. It is, however, not 
without problems. In particular, the predicted galaxy stellar 
mass function exhibits an excess of low-to-intermediate mass 
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galaxies at all redshifts, the fraction of red galaxies among 
low mass galaxies is too high compared to observational 
data, and the mo del over-predicts the clustering signal of 
faint galaxies (e.g. Weinmann et ahlliood : IWang et"aLll2008l : 
iFontanot et all |2009|) ." It is important to note that these 
'failures' are not specific of the particular model adopted 
in this study. Rather, they appear to be common to most 
(all) models that have been published recently, and have 
likely related causes. No satisfactory explanation and/or so- 
lution has been found yet to solve all of th e problems men- 
tion ed above (for recent attempts, see e.g. IGuo et al]|201ll 
and IWang. Weinmann fc Neisteinll201ll ') . 

The model used in this study neglects environmental 
physical processes such as ram-pressure and harassment, 
but assumes that when galaxies are accreted onto a more 
massive system, the associated hot gas reservoir is stripped 
instantaneously. This induces a very rapid decline of the 
star formation histories of satellite galaxies, and contributes 
to create an excess of red a nd passive galaxie s with re- 
spect to the observations (e.g. IWang et al.ll2007l V In recent 
studies, a more gr adual stripping of the hot gas reservoir 
has been assumed (iKang fc van den Boschll2008l : [ptant et al.l 
l2008l : IWeinmann et al.ll2010l : IGuo et al.ll201ll). following re- 
sults fro m numerical simu lations ( McCarthv et al.ll2008l . but 
see also ISaro et al.|[2010l ). Albeit improved, the agreement 
with observational measurements is far fro m satisfactory 
(e.g. IGuo et al.ll201ll : [Weinmann et al.ll201ll '). In this study, 
we will make limited use of the physical properties of galax- 
ies as predicted by the model. In particular, we will use 
only the predicted galaxy stellar mass, and we will focus 
on the resulting galaxy merger trees and their dependence 
on the parent halo mass. In the following, we will focus on 
haloes selected from the box of the simulation correspond- 
ing to z=0, with mass larger than ~ 10^'^ M©. The halo 
mass is computed from the A^-body simulation, as the mass 
within a sphere enclosing a mean overdensity that is equal 
to 200 times the critical density of the Universe (M2oo). 
At the resolution of the Millennium Simulation, haloes with 
M200 ~ 10^'^ M0 contain on average about 15 galaxies within 
the virial radius (-R200 in this study), with stellar mass larger 
than > 10' M0 . To investigate trends as a function of halo 
mass, we have randomly selected fifteen haloes in each of 
the following mass bins: ~ 10^^, ~ 5 x 10^^, 10^*, and 
~ 5 X 10^^ M0. These haloes contain 9201 satellite galax- 
ies within their virial radii, that represent the basic sample 
of our analysis. 6929 of these satellites have present day stel- 
lar mass between lO" and 10^° Mq, 2081 have mass between 
10^" and 10^^ M!q, and only 191 of them have stellar mass 
larger than lO" Mq. 

In order to clarify the method and the definitions that 
we have adopted and that we will use in the following, we 
show in Figure [T] an example of a galaxy merger tree. The 
galaxy considered is shown at the top of the plot, and all its 
progenitors (and their histories) are shown going backward 
in time. The size of the symbols scales with the galaxy stel- 
lar mass, while different colours are used for central galaxies 
(shown in black), satellite galaxies associated with a dis- 
tinct dark matter substructure (green), and satellite galax- 
ies whose parent dark matter subhaloes have been stripped 
below the resolution of the simulation (these are shown in 
red). We note that only central galax;ies are fuelled by cool- 
ing flows in our model. The leftmost branch in Figure [T] is 



" - M..„ = 7.90 X 10'° M„ 

- M|,„„ = 5.53 X lO'^^M, 
^Mv = -21.1130 

2 _ B - V = 0.913095 

- type - 1 




Figure 1. Example merger tree for a galaxy of stellar mass 
~ 8 X 10^° Mq, residing in a halo of mass ~ 5.5 X 10^* Mq. 
The size of the symbols scales with the galaxy stellar mass while 
different colours are used for central galaxies (black), satellite 
galaxies associated with distinct dark matter subhaloes (green), 
and satellites whose parent substructures were stripped below the 
resolution of the simulation (red). The empty black circle at look- 
back time ~ 5.5 Gyr corresponds to the lookback time before the 
galaxy is accreted onto the main progenitor of the final halo. 



obtained by connecting the galaxy to its most massive pro- 
genitor (often referred to as the mam progenitor) , at each 
node of the tree. 

In the following, when tracing the history of our model 
galaxies, we will effectively refer to the main progenitor of 
the galaxy at each time. In addition, we define two different 
characteristic times: (i) the lookback time before a galaxy 
becomes satellite of a larger halo, and (ii) the lookback time 
before the galaxy is accreted onto the main progenitor of 
the final halo. In the following, we will often refer to the for- 
mer characteristic time as tsa,t or the time when the galaxy 
becomes a satellite, and to the latter as thaio or the time 
when the galaxy is accreted onto the final group/cluster. In 
the example shown in Figure [TJ tsat corresponds to the last 
time when the colour of the main progenitor is black be- 
fore becoming green (~ 7 Gyr), while thaio is marked by 
an empty black symbol (~ 5.5 Gyr for the example con- 
sidered). As we will discuss in detail in the following, the 
adoption of different characteristic times is l argely respon- 
sible f or t he different conclusio ns reached by iBerrier et al.l 
(| 2009l ) and lMcGee et al.l (|2009l ). 

By construction, the difference between tsat and thaio 
has to be positive, or zero. Figure [5] shows the distribution 
of tsat — ihaio, for galaxles of different stellar mass. As the 
trends are similar for haloes of different mass, we have con- 
sidered here all galaxies in all haloes analysed for this study. 
The figure shows that for ~ 40 to ~ 60 per cent of the galax- 
ies, the two events considered happen at the same time. For 
the remaining galaxies, these two events can differ by up to 
~ 8 Gyr. A larger fraction of the most massive galaxies in 
our sample are accreted as centrals and, for these galaxies, 
the difference between the two events considered is on aver- 
age smaller than for less massive galaxies. This results from 
the combination of different factors. There is a strong cor- 
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Figure 2. Distribution of tlie difference between the lookback 
time corresponding to the event of being accreted onto another 
system (i.e. becoming a satellite), and that corresponding to the 
accretion onto the final group/cluster (see text for details). Differ- 
ent colours correspond to galaxies of different present day stellar 
mass, as indicated in the legend. 



If one defines as clusters all haloes more massive than 
~ 10" M0, the fi gure shows that only about 10 per cent of 
the cosmic galaxy population resides in clusters at present, 
and that this is approximately independent of the galaxy 
stellar mass threshold considered. As expected, this fraction 
decreases at higher redshift, as massive clusters form rela- 
tively late. The figure also shows that more massive galaxies 
tend to reside in more massive haloes: about 50 per cent of 
the galaxies with stellar mass > 10^^ M0 reside in haloes 
more massive than > 10^^ M©, that are sometimes referred 
to as groups in the literatur^l. When considering all galaxies 
with stellar mass > 10^" Mq, however, this fraction drops 
to only about 30 per cent at the present day. It is even lower 
for lower stellar mass limits. Therefore, it is not generally 
true that 'most galaxies reside in groups', if by groups one 
refers to haloes with mass > lO'^'^ M©. As just discussed, this 
statement depends significantly on the stellar mass limit of 
the observed sample. We stress that the definitions adopted 
in this paper are theoretical ones, and that they do not nec- 
essarily always represent a good proxy for observational def- 
initions. In this paper, we will work in this 'theory space', 
and use the mass of the parent halo (specifically, M200) as 
a proxy for the environment. In future work, we plan to 
extend our analysis beyond the virial radius, and to use en- 
vironmental definitions that are closer to those commonly 
adopted in the literature. 



relation between the stellar mass of central galaxies and the 
parent halo mass, and massive haloes form relatively late in 
hierarchical cosmogonies. So the most massive galaxies can 
only be satellites of very massive haloes: the main progen- 
itor of the parent halo has to be already quite massive in 
order to accrete a second system that is also large enough 
to contain a massive galaxy. In addition, massive galaxies 
that were accreted early on during the history of the parent 
halo, will be rapidly dragged closer to the inner regions by 
dynamical friction, and eventually merge with the central 
galaxy of the final halo. Therefore, these galaxies will not 
appear in our sample that only contains satellite galaxies 
surviving at redshift zero. 

Before studying the history of the environments of our 
model galaxies, it is interesting to analyse their present day 
distribution as a function of the parent halo mass, and how it 
varies as a function of cosmic time. FigureOshows the differ- 
ential (top panels) and cumulative (bottom panel) distribu- 
tion of galaxies for three different galaxy stellar mass thresh- 
olds (increasing from left to right columns). These have been 
obtained by using the entire volume of the Millennium Sim- 
ulation, and all galaxies with stellar mass > 10® Mq. Differ- 
ent linestyles correspond to different cosmic epochs (solid is 
for present day, dotted for 2 ~ 0.5, and dashed for 2: ~ 1). 
Each distribution has been normalized to the total number 
of galaxies above the stellar mass limit indicated in the leg- 
end. Note that the cumulative distributions in the bottom 
left and middle panels do not reach unity at the lowest halo 
mass limit shown: ~ 55 per cent of the galaxies with stel- 
lar mass larger > 10® Mq reside today in haloes with mass 
< lO^'^ Mq. For galajcies more massive than ~ 10^" M©, the 
corresponding fraction drops to ~ 40 per cent, and it be- 
comes zero for galaxies more massive than ~ lO'^^ Mq. 



3 THE HISTORY OF GROUP AND CLUSTER 
GALAXIES 

In order to study the environmental history of our model 
satellite galaxies, we have traced back their histories and 
stored the values of their parent halo mass at all times. In 
addition, for each galaxy, we have computed the time when 
it becomes a satellite, and the time when it is accreted onto 
the main progenitor of the final halo, and stored the values 
of the parent halo mass corresponding at these times. In the 
following, we will refer to the value of the parent halo mass 
at tsat as Mhaio(sat), and to the corresponding value at thaio 
as Mhaio(halo). 

Figure |4] shows the projected comoving trajectories of 
three cluster galaxies, re-centred with respect to the posi- 
tions of the main progenitor of the final halo at each time. 
The galaxies have been selected randomly between the satel- 
lites of haloes with mass ~ 10^* M©, and with present stel- 
lar mass > 10" Mq. The progenitors of these galaxies come 
from extended regions around the main progenitor of the 
final halo, and their parent halo mass (encoded in the size 
of the symbols) increases progressively as the galaxies get 
closer. Figure [5] shows the mean evolution of the stellar mass 
(dashed lines) and parent halo mass (solid lines), computed 
considering all satellites of 15 haloes with mass ~ 10" Mq 
(qualitatively, these results do not change significantly when 
considering different ranges of parent halo mass) . The mean 
evolution of the galaxy stellar mass is very similar in the 
three mass bins considered. The evolution of the parent halo 

Note t hat the Local Gro up, with an estimated mass of 5 X 
10^^ Mq llLi fc White! I2OO8I I, would not be classified as a group 
according to this definition. 
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Figure 3. Differential (top panels) and cumulative (bottom panels) distributions of galaxies as a function of their parent halo mass. 
Different columns correspond to different galaxy stellar mass thresholds, as indicated in the legend. Different linestyles correspond to 
different cosmic epochs (solid, dotted and dashed lines are used for z ~ 0, 0.5, and 1, respectively). Each distribution is normalized to 
the total number of galaxies above the stellar mass limit considered (without dividing by the bin size). 
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Figure 4. Projected (comoving) positions of the main progenitors 
of three cluster galaxies, selected randomly between the satellites 
with stellar mass ^ 10^^ Mq, residing in haloes with mass ~ 
10^* Mq. Coordinates are re-centred with respect to the positions 
of the main progenitor of the final cluster. The size of the symbols 
scales with the parent halo mass, while the colours encode the 
lookback time. 
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Figure 5. Mean evolution of the galaxy stellar mass (dashed 
lines), and of the parent halo mass (solid lines). These have been 
computed by averaging the evolution of all satellites of haloes 
with mass ~ 10^^ Mq. Different colours correspond to galaxies of 
different present day stellar mass, as indicated in the legend. 



masses (solid lines) shows clearly that more massive galax- 
ies are sitting on average in more massive haloes for most of 
their evolution (before starting being accreted onto the final 
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haloes). The difference in mass is, however, not large. In ad- 
dition. Figure [5] also shows that the most massive galaxies 
surviving at redshift zero fell onto the main progenitor of 
the final cluster later than their less massive counterparts: 
all lines superimpose when gala^xies are in the same halo, and 
for the most massive galaxies this happens only at z ~ 1. 

It is very instructive to look at some properties of 
group/cluster satellites at the characteristic times defined 
above. The left panel of Figure [6] shows the time since the 
galaxy became a satellite as a function of the corresponding 
parent halo mass, for all galaxies residing in a single cluster 
of mass M 5 X 10^* M0 at ^ = 0. The solid line shows 
the mass accretion history of the parent cluster, while the 
size of the symbols scales with the present day galaxy stel- 
lar mass, as indicated in the legend in the right panel. The 
figure shows that cluster members with low stellar mass be- 
come satellites (are accreted onto more massive haloes) at 
all times, while the (few) most massive galaxies that sur- 
vive as satellites at redshift zero tend to become satellites 
later compared to the overall cluster galaxy population. As 
discussed in the previous section, this is expected in a hier- 
archical cosmology. 

Since the number of galaxies in a halo increases with 
halo mass, those galaxies that became satellites when sit- 
ting in massive systems also have larger numbers of com- 
panion galaxies at that time, as shown in the right panel 
of Figure [6] Interestingly, the bottom right corners of both 
panels in Figure |6] appear to be rather 'empty'. The mass 
of the cluster considered in our particular example grows 
significantly down to a lookback times ~ 4 Gyr (that corre- 
sponds to 2: ~ 0.4), through the accretion of few relatively 
massive structures. At lower redshift, the mass growth of 
the parent halo is not significant, and it is mainly driven by 
the accretion of small haloes and diffuse material. Most of 
the galaxies that become satellites at lookback times lower 
than ~ 4 Gyr are accreted directly onto the main progen- 
itor of the final cluster, and have low stellar mass. As we 
will see below, these results hold for the statistical sample 
considered in this study. It is interesting that this regime 
of 'slow mass growth' for the parent haloes coincides with 
the cosmic epoch that witnesses the most striking transfor- 
mations for the cluster galaxy population, both in terms of 
their star form ation activi ty and morpho logical mix (e.g. 
iButcher fc Oemle r 1984; Desai et al.ll2007l '). 

Figure [7] shows the distributions of the times when 
group/cluster members become satellite galaxies (top left 
panel), of the times when the galaxies are accreted onto 
the final halo (bottom left panel), and the corresponding 
distributions of parent halo mass (right panels). As there 
is no significant trend as a function of final halo mass, we 
have stacked together all haloes in our sample. The top left 
panel of Figure[7]shows that cluster/group members became 
satellite galaxies over a wide range of lookback times. The 
distributions for low and intermediate mass galaxies exhibit 
a peak at very early epochs (~ 10 Gyr), while the distri- 
bution obtained for the most massive galaxies is peaked at 
later times, with most of them becoming satellite between 4 
and 9 Gyr ago. 

The top right panel of Figure [7] shows a strong corre- 
lation (albeit with a relatively large scatter) between the 
present day galaxy stellar mass and the parent halo mass at 
iaat- This refiects the strong correlation between halo mass 
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Figure 8. Median (solid lines) and percentiles (25th and 75th) 
of the lookback times since cluster galaxies became satellites, as 
a function of their present day (three-dimensional) distance from 
the cluster centre. This has been obtained by stacking all haloes 
with mass ~ 5 X 10^^ Mq in our sample. Different colours refer 
to different galaxy stellar mass bins, as indicated in the legend. 



and galaxy stellar mass for central galaxies. This correla- 
tion is preserved for satellite galaxies because their stellar 
mass does not increase significantly after being accreted, as 
their star formation is efficiently suppressed over a quite 
short time-scale in our model. The distribution of parent 
halo masses found for the most massive galaxies is also the 
widest, due to the fact that the relation between the galaxy 
stellar mass and the parent halo mass flattens at high stellar 
masses. The distributions shown in the top panels of Figure[7] 
look different when one considers the time when the galaxies 
were accreted onto the most massive progenitor of the final 
halo (thaio)- These distributions are shown in the bottom 
panels of Figure [7] Since, as discussed above, the two times 
considered tend to coincide for the most massive galaxies, 
the distributions computed for these galaxies are not signif- 
icantly different. Results are instead different for low and 
intermediate stellar mass galaxies. A large fraction of these 
tend to be accreted when they are already satellite galaxies, 
so that the distributions of thaio are shifted towards later cos- 
mic epochs with respect to the corresponding distributions 
of tsat- As a consequence, if one looks at the distribution of 
halo masses at thaio, this is found to be 'bimodal' for these 
galaxies. The peaks at lower halo masses are made up of 
galaxies that are accreted when they are central galaxies of 
their own haloes, while the peaks corresponding to larger 
parent halo masses are populated by galaxies that are ac- 
creted as satellites of relatively large systems. In particular, 
we find that ~ 48 per cent of the galaxies in the lowest stel- 
lar mass bin considered were accreted as satellite galaxies. 
The fraction is slightly lower (~ 43 per cent) for the inter- 
mediate mass bin considered, while it drops to only ~ 23 
per cent for the most massive galaxies in our sample. 

Substructures that were accreted onto the main progen- 
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Figure 6. Left panel: lookback time corresponding to the event of becoming a satellite galaxy, as a function of the parent halo mass at 
the same time. Each symbol corresponds to a galaxy, with the symbol size scaling proportionally to the present day stellar mass. The 
example shown corresponds to a single halo of mass Af ^ 5 X 10^* Mq, selected at 2 = 0. The solid thick line shows the mass evolution 
of the cluster main progenitor. Right: lookback time corresponding to the event of becoming a satellite as a function of the galaxy's 
number of satellites at the same time. 



itor of their parent halo at early times had relatively short 
orbital periods. So these should be located, on average, in the 
inner regions of the final halo. In addition, these haloes will 
have suffered from dynamical friction for a longer time with 
respect to haloes of similar mass but accreted later. These 
two factors combine to create a strong correlation between 
the accretion time an d the cluster-cen tric distance of dark 
matter substructures (|Gao et al.ll2004l ). Figure [8] shows that 
this correlation is also found for cluster galaxies, although 
with a large scatter. Therefore, a radial dependence of galaxy 
properties is, at least in part, a natural consequence of the 
fact that mixing of the galaxy population is incomplete dur- 
ing cluster assembly. The weakest trend is found for the 
most massive satellites surviving at redshift zero that, as 
discussed above, were accreted later than their less massive 
counterparts. The radial trend is stronger for the low and 
intermediate stellar mass galaxies considered: those that are 
located closer to the centre have been orbiting as satellite 
gala:xies for about 11 Gyr, while those that are located at the 
outskirts of the haloes have been satellites for about 6 Gyr 
on average. As we will discuss in more detail in the following 
section, the trends shown in Figure [H] can be combined with 
the observed fractions of red/passive galaxies as a function 
of cluster-centric distance to put constraints on the physical 
processes (and corresponding timescales) responsible for the 
suppression of star formation activity in group and cluster 
galaxies. 

A figure similar to our Fig ure [8] has been shown in 
IWeinmann. van den Bosch fc Pas guali (2011) who took ad- 
vantage of the same simulated galaxy catalogues but used 
a cluster sample selected mimicking a specific observational 
data-set, and plotted results as a function of the projected 
cluster-centric radius. Simil ar findi ngs have also been dis- 
cussed in a recent work by iSmith et al.l (|201^ ') . These au- 



thors have based their analysis on publicly ava i lable cata- 
logues from the model presented in 'Font et al.' ('2008', also 
based on the Millennium Simulation). In particular, they 
have selected all galaxies residing in the five most massive 
haloes in the simulation, and studied how various character- 
istic epochs of the evolution of model galaxies relate to their 
present da y projected d i stanc e from the cluster centre. The 
results by ISmith et al.l ||2012|') are in very good agreement 
with those by Weinmann et al.l (|201ll ) and with ours, con- 
firming that these findings are not significantly dependent 
on the particular model used, but mainly driven by structure 
formation. 



4 HISTORY BIAS 

The results shown in the previous section demonstrate that 
the galaxy population of groups and clusters is character- 
ized by a variety of accretion histories, and that these are 
strongly dependent on the final galaxy stellar mass. In this 
section, we will discuss how the results found above can in- 
fiuence the observed environmental trends. In particular, we 
will focus on two observables: the fraction of passive/red 
galaxies that reside in haloes of different mass, and its vari- 
ation as a function of the cluster-centric distance. 

Figure[5]shows the fraction of galaxies that spent a time 
longer than that given on the x-axis in haloes with mass 
larger than 10^^ (cyan), lO" (green), and 10^* M© (orange). 
In this Figure, each solid thin line corresponds to a different 
halo, while the thick solid lines show the mean obtained for 
the 15 haloes considered in each sample. As a consequence of 
the hierarchical structure formation, the fraction of galaxies 
that have resided in haloes more massive than Mhaio for a 
time longer than Thaio decreases when considering both more 
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Figure 7. Distributions of the times when group/cluster members become satellite galaxies (top left panel), and of the times when the 
galaxies are accreted onto the main progenitor of the final halo (bottom left panel). The right panels show the corresponding distributions 
of parent halo masses. The distributions shown have been computed using all haloes selected at 2 = 0. Different colours correspond to 
different present day galaxy stellar mass, as indicated in the legend. 



massive haloes and longer times. So it is not surprising that 
the thick green lines are always below the thick blue lines, 
and that these are always above the red lines. Interestingly, 
the figure shows a relatively large halo-to-halo scatter: e.g. 
focusing on the most massive haloes in our sample, the frac- 
tion of galaxies that have spent more than 6 Gyr in haloes 
more massive than 10" Mq varies between ~ 20 per cent 
and ~ 80 per cent. If the time spent in a more massive 
halo can be related to the star formation activity (or lack 
thereof) of a galaxy, this halo-to-halo scatter can be linked 
to the different fractions of passive galaxies measured in dif- 
ferent clusters (e.g. in the nearby Coma and Virgo clusters, 
IWeinmann et al1l201ll ) and, more in general, to the scatter 
measured for some properti es of the galaxy populations of 
group s and clusters (see e.g. IPe Lucia. Fontanot fc WilmanI 
|2012| ). Interestingly, the observed fraction of red/passive 
galaxies does indeed show a large halo-to-h alo scatter, that 
appears to increase at lower halo masses (|Poggianti et al.l 
120061 : iBalogh fc McGeell2010l ). It should be noted, however, 
that observational uncertainties are usually much larger for 



lower mass systems. In addition, as we will show in the 
next sections, the fraction of passive galaxies increases only 
weakly as a function of halo mass, and is larger than ~ 0.4 in 
haloes of mass ~ 10^'^ M©. Taken at face value, these results 
favour a value for A/haio (i.e. the halo mass above which star 
formation stars to decline) lower than 10^* M0 that would 
produce a too low fraction of passive galaxies in haloes with 
present day mass ~ 5 x 10^"^ Mq or lower. 

We discuss this in more detail in the following sub- 
sections, where we investigate how the trends shown in Fig- 
ure [9] depend on galaxy stellar mass and on the cluster- 
centric distance, and how the theoretical expectations 
compare to the estimated fractions of 'passive' or 'red' 
group/cluster galaxies. 



4.1 Observational data 

To constrain our theoretical models, we take advantage 
of the group catalogue by lYang et al.l ||2007]) comple- 
mented with the cluster catalogue by von der Linden et al.l 
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Figure 9. Fraction of galaxies that spent a time longer than that given on the x-axis in haloes with mass larger than 10^^ (blue), lO'^^ 
(green), and lO^"* Mq ( orange). Each solid thin line corresponds to a single halo, while the corresponding thick lines represent the mean 
of the distributions. Different panels correspond to galaxies residing in haloes of different mass today (increasing from left to right), as 
indicated in the legend. 



both based on DR4. Specifically, we use the sam- 
ple II of the Yang e t al. g roup catalogue, as described in 
Ivan den Bosch et ahl (|2008l n. We use only satellite galax- 
ies (as in the theoretical predictions discussed below), and 
weight results according to the maximum value out to which 
they can be observed, to account for Malmquist bias. Our 
sample of satellite galaxies based on the Yang et al. cata- 
logue contains 42363 galaxies residing in haloes w ith mass 
larger than 10^^ Mq . As for the von der Linden et al.i (i2007l ) 
catalogue, their final sample consists of 625 systems at red- 
shifts between 0.03 and 0.1, with masses between ~ 10^^ and 
~ 10^^ Mq. In the following, we will consider only the 214 
systems more massive than ~ lO^'' Mq from this catalogue. 

In addition, we make use of the data catalogues for 
SDSS DR7 from MPA/JHU to obtain up-to-date estimates 
for stellar masses and specific star formation rates (SS- 
FRsfl. Stellar masses are estimated using fits to the pho- 
tometry, a nd are statistically in ag reement with the esti- 
mates from iKauffmann et al.1 hooj). The SSFR s are based 
on those published in Brinchmann et al] (|2004l ). but with 
several modifications regarding the treatment of dust atten- 
uation and aperture corrections. AH details can be found at 
the reference w ebpage. 

F ollowing IWeinmann et al.l (|2006l ) and iKimm et al.l 
(|2009l ). we have adopted the following demarcation line be- 
tween red and blue galaxies: 



\9 



0.7 



0.032 C'^Mr 



51og/i-f 16.5) 



(1) 



All galaxies with l og(SSFR) < —1 1 are c lassified as pas- 
sive. As shown bv IWeinmann et al] l|2010l ). this cut corre- 
sponds roughly to the location of the minimum in the bi- 
modal distribution of SSFRs in the observations. For com- 
parison, we have also used alternative estimates of the SS- 



^ The group catalogue is publicly available 
|http://www.astro.umass.edu /~xhyang/ G roup, html | 

" publicly available 



The 



data 



at: 



at: 



|http://www.mpa-gar ching.mpg.de/SDSS71 



FRs base d on UV data from the GALEX satellite, as de- 
scribed in iMcGee et al.l (|201ll ). This results in a reduction 
of the sample size by a factor of about 10, as not all SDSS 
galaxies have GALEX coverage. The scatter in the obser- 
vational measurements given below gets larger, but quali- 
tatively results do not change. Therefore, we only show our 
measurements of passive fractions based on the Brinchmann 
et al. estimates of the SSFRs. 

Finally, the cluster sample is complemented with the 
d atasets for Coma and V irgo that are described in detail 
in lWeinmann et al.l (|201ll '). For these datasets, stellar mass 
estimates are obtained from fits to the photometry, using the 
S DSS g- and r-band filters , acco rding to the fitting formula 
bv lZibetti. Chariot fc Rb] (|2009l ). 



4.2 Cluster Radial Dependence 

Figure [TO] shows the fraction of model galaxies that spent 
more than 3, 5, 7, and 8 Gyr (lines of different colour) in 
haloes more massive than the thresholds indicated in the 
legend of each panel (increasing from left to right). Top 
and bottom panels are for galaxies with stellar masses in 
the ranges log(Mstar/M0) = [9 - 10] and log(Mstar /Mq) = 
[10 — 11], respectively. Theoretical predictions have been 
obtained averaging all haloes in our sample with mass 
larger than 10^* Mq, and projecting cluster members on 
the xy plane. Red diamonds, blue triangles, and black cir- 
cles show the fractions of red galaxies for the Coma clus- 
ter, for the Virgo clu s ter, a nd for the cluster catalogue by 
Ivon der Linden et al.l l|2007l ). respectively. Finally, red cir- 
cles show the fraction of passive galaxies estimated using 
the catalogue by von der Linden et al. as in lWeinmann et aU 
(2010). Error bars denote confidence intervals corresponding 
to a probability of 68 per cent (la) from qu artiles of the bet a 
distribution (for details on the method, see ICamero nl l201ll '). 

Interestingly, the r ed fractions estimated from t he clus- 
ter sample defined in Ivon der Linden et al.l (|2007l ) are in 
quite good agreement with the red fractions estimated for 
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Figure 10. Fraction of model galaxies that spent more than 3, 5, 7, and 8 Gyr (lines of different colour) in haloes more massive 
than 10^^ (left panel), 10^^ (middle panel), and 10^* (right panel). Top panels are for galaxies with stellar mass in the range 
log(Matar/M0) = [9 — 10], while bottom panels are for galaxies in the range log(Matar/MQ) = [10 — 11]. Theoretical predictions have 
been obtained averaging all haloes in our sample with mass larger than lO'^* Mq , and projecting cl uster members on the xy plane. Red 
diamonds and blue triangles show the red fractions for the Coma and Virgo clusters estimated bvlW einmann et al.l ll201lh . Black and 
red circles show the red and passive fractions estimated using the cluster catalogue bv lvon der Lindeii'et'^Lri|200^ TrresDectivelv. Error 
bars denote confidence interval estimates measured as described in the text. 



the Coma and Virgo clusters. This suggests that neither 
contamination by background galaxies (this affects at some 
level the Coma cluster catalogue as it includes photomet- 
ric members), nor incompleteness due to fiber completeness 
(which affects the cluster catalogue by von der Linden et 
al.) influence our results significantly, at least out to a pro- 
jected radius of ~ 1 Mpc. Both the observed red and passive 
fractions exhibit a clear dependency on stellar mass: more 
massive galaxies exhibit a shallower radial trend with re- 
spect to their less massive counterparts (compare bottom 
and top panels of Figure [10)) . in qualitative agreement with 
the trends shown in Figure [8] The measured fractions of 
passive galaxies are generally lower than the fractions of red 
galaxies, indicating that a fraction of all star forming galax- 
ies are red because dusty. The offset between the red and 
passive fractions depends on galaxy stellar mass and (albeit 
weakly) on the projected distance from the cluster centre. 
In particular, using the res ults based on the catalogue by 
Ivon der Linden et al.l l|2007t l. we find that a fraction vary- 



ing between ~ 5 per cent (at small projected distances) 
to ~ 10 per cent (at the largest projected distances con- 
sidered) of the red galaxies are star forming dusty systems 
in the stellar mass bin log(Mstar/Af0) = [10 — 11]. For the 
lower mass bin considered, the fraction of star forming dusty 
galaxies gets larger, varying between ~ 20 per cent close to 
the cluster centre, to ~ 24 per cent at projecte d distances 
of ~ 1 Mpc. Our results contradict those by IWolf et al.l 
(|2009h who analyse the properties of dusty red galaxies in 
the A901/2 cluster complex at 2: ~ 0.17, and argue that 
dusty star forming galaxies are rare for log(Mstar/Af0) < 10 
while appearing predominantly in the stellar mass range of 
log(Afstar/M0) = [10 - 11]. We note, however, that Woff et 
al. defined their population of dusty red galaxies differently 
than done in our study. In particular. Wolf et al. defined 
red dusty galaxies from a full spectral energy distribution 
fit to the medium-band photometry of the COMBO-17 sur- 
vey. This will be more uncertain for faint galaxies (i.e. for 



© 2011 RAS, MNRAS OOO.ITlfTTl 



The environmental history of group and cluster galaxies 11 



low mass galaxies), a caveat that applies to the passivity 
criterion used in our study as well. 

The theoretical predictions shown in Figure[lO]are shal- 
lower when considering lower values of Thaio at fixed Mhaio 
values, and for the largest values of Mhaio considered (right 
panels). Again, this is a consequence of hierarchical struc- 
ture formation: more massive haloes are formed later than 
their less massive counterparts so that galaxies could spend 
less time in relatively massive systems. A comparison be- 
tween theoretical predictions and observational data shows 
that different combinations of Thaio and Afhaio can provide 
a relatively good agreement with the observed radial trends 
of passive galaxies. For example, for galaxies with stellar 
mass in the range log(Mstar/MQ) = [9 — 10] (upper pan- 
els), the observed trend is very close to that obtained when 
considering galaxies that have been sitting in haloes more 
massive than 10^'^ Mq for more than ~ 8 Gyr. A quali- 
tatively good agreement (but with a slope shallower than 
observed) is obtained when considering galaxies that have 
been sitting in haloes more massive than 10^"* M© for more 
than ~ 5 — 7 Gyr. For the more massive bin considered 
(log(Matar/M0) = [10— 11]), the fraction of passive galaxies 
has a significantly shallower radial trend. Qualitatively, it 
is very similar to the trend obtained considering the frac- 
tion of galaxies that spent more than ~ 7 Gyr in haloes 
more massive than 10^^ M0, but also to the trend obtained 
considering galaxies that spent more than ~ 5 Gyr in haloes 
more massive than 10^^ M0, or more than ~ 3 Gyr in haloes 
more massive than lO^'' . 

4.3 Halo mass dependence 

Another important and independent constraint to our theo- 
retical predictions is given by the observed fraction of red or 
passive galaxies as a function of the present day halo mass. 
Figure [TT] compares the fraction of model galaxies that spent 
more than 3, 5, and 7 Gyr (from left to right panels) in haloes 
more massive than 10^^ 10^^ and lO" Mq (lines of different 
style) with observational estimates of red and passive galax- 
ies (top and bottom panels, respectively). To construct the 
datasets sho wn in this fi g ure, w e have used the DR4 group 
catalogue bv lYang et al.l l|2007l ). as detailed above. 

Figure [TT] shows that both the fraction of red galax- 
ies (top panels) and that of passive galaxies (bottom pan- 
els) tend to increase as a function of halo mass. In addi- 
tion, at fixed halo mass, the fraction of red/passive galax- 
ies is larger for more massive galaxies. The res ults obtained 
are in good agreement with those presented in lKimm et al.l 
|200a ). As noted above, the fraction of passive galaxies is 
generally lower than that of red galaxies. In particular, we 
find that about 16 per cent of the galaxies are red but 
star forming for the lowest galaxy stellar mass bin con- 
sidered (log(Mstar/A/Q) ~ 9.75). The fraction reduces to 
about 5 per cent for the largest stellar mass bin considered 
(log(Mstar/M0) ~ 11.25): i.e. a larger fraction of red mas- 
sive galaxies are truly passive, while for less massive galax- 
ies the fraction of star forming but dusty (and therefore 
red) objects increases. This is in agreement with what dis- 
cussed about Fig ure IIOI and in contrast with results from 
IWolf et all (|2009l i. 

In our models, more than 90 per cent of the galaxies 
with stellar mass larger than log(Mstar/M0) ~ 9.75 have 



spent more than 3 Gyr in haloes more massive than 10 M© 
(dotted lines in the left panels of Figure lllfl . When consid- 
ering longer times (middle and right panels) the fractions 
decrease, and one can see some dependency on the galaxy 
stellar mass (larger fractions are obtained for more mas- 
sive galaxies), in particular for the longest times considered. 
The fractions of galaxies that have spent equal time in more 
massive haloes are lower, as a natural consequence of hier- 
archical structure formation. In particular, our model pre- 
dicts that only about 60 per cent of the galaxies in the most 
massive haloes considered in our sample have spent more 
than 7 Gyr in haloes more massive than lO'^^ Mq (dashed 
lines in right panels). The fraction decreases rapidly for less 
massive haloes. For the most massive galaxies considered 
(red symbols), the observed passive fractions are in quite 
nice agreement with the estimated fractions of galaxies that 
spent more than 5 Gyr in haloes more massive than 10^'^ M©. 
For less massive galaxies, longer times seem to work better 
(compare cyan lines and points in the bottom middle and 
right panels of Figure llip . 

4.4 Central quenching and Satellite quenching 

Putting together both the estimated fractions of passive 
galaxies as a function of halo mass and those as a func- 
tion of cluster-centric distance, there seems to be a rough 
timescale/halo mass combination that gives results in quite 
good agreement with observational measurements: Mhaio ~ 
10^'^ Mq and Thaio ~ 5 — 7 Gyr is needed for a galax;y to 
be quenched with close to 100 per cent probability, with a 
tendency for slightly longer timescales and slightly larger 
halo masses for lower stellar mass galaxies. However, the 
interpretation of the last two figures shown is not trivial 
because one has to consider that there is a continuous evo- 
lution of the halo mass hosting the galaxies that end up in 
a group/cluster. In addition, it is necessary to account for 
galaxies that become passive because of 'internal' physical 
processes {central- quenching), rather than for 'environmen- 
tal' processes {satellite quenching). 

The importance of central-quenching can be estimated 
by consideri ng the fractio n of re d /passive central galaxies. 
Figure 3 by iKimm et al.l (|2009l ) shows how this fraction 
varies as a function of the galaxy stellar mass. Virtually all 
central galaxies in the first stellar mass bin considered in our 
study, log(Mstar/Af0) = [9.5 — 10], are active. For the larger 
stellar mass bins that we have considered in Figure [TT] the 
fraction of passive central galaxies increases to ~ 10 — 20 
per cent for the bin log(Mstar/MQ) = [10 - 10.5], and to 
~ 40 - 50 per cent for the bin log(Mstar/M0) = [10.5 - 11]. 
For the most massive galaxies considered, ~ 70 per cent of 
the central galaxies are already passive. We can therefore 
distinguish two galaxy stellar mass regimes: 

(i) an satellite quenching dominated regime: galaxies with 
log(Afstar/M0) < 10.5 had their star formation rate mainly 
suppressed by the environment. Our results show that the 
estimated fractions of passive galaxies in this mass regime 
are in quite good agreement with the expected fraction of 
galaxies that have spent more than ~ 5 — 7 Gyr in haloes 
more massive than 10^'' M0 ; 

(ii) a central quenching dominated regime: galaxies with 
log(Mstar/M0) > 10.5 whose star formation rate has been 
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Figure 11. Fraction of model galaxies tiiat spent more than 3 (left panel), 5 (middle panel), and 7 Gyr (right panel) in haloes more 
massive than 10^^ (dotted lines), 10^'^ (dashed lines), and 10^* Mq (solid lines). Fractions are plotted as a function of the present day 
halo mass. Different colours correspond to different galaxy stellar mass bins: cyan corresponds to log(Mstar/M0) = [9.5 — 10], green to 
log(Mstar/M0) = [10-10.5], orange to log(Mstar/M0) = [10.5 - 11], and red to log(Mstar/M0) = [11 - 12]. Data points wit h error bars 
are ob servational estimates of red fractions (top panels) and passive fractions (bottom panels) based on the DR4 catalogue by I Yang et al.l 
Error bars denote confidence intervals from quartiles of the beta distribution. 



reduced primarily by internal physical processes. For these 
galaxies, the estimated passive fractions shown in Figure [TT1 
favour a somewhat lower value of Thaio ~ 5 Gyr. This could 
be due to the fact that these galaxies have spent some frac- 
tion of their lifetime as central galaxies. 

In order to determine the importance of environmental 
processes in suppre ssing the star formation rate s of satellite 
galaxies, we follow Ivan den Bosch et al.l (|2008l ) and define 
the following transition fraction: 

f fraC red, sat fraC rcd.ccn / 

traCcnv = 7 ■ — 

iraC bluc.ccn 

where frac red, sat is the fraction of satellite galaxies that are 
red, frac red.cen is the fraction of central galaxies that are red, 
and frac biue,cen is th e fraction of central galaxi es that are 
blue. As discussed bv lvan den Bosch et al.l (|2008t ). under the 
assumption that the present day population of central galax- 
ies is representative of the progenitors of present-day satel- 
lite galaxies of the same stellar mass, the numerator of Eq.[2] 



gives the fraction of all satellite galaxies that have undergone 
a blue-to-red transition after their accretion. Dividing by the 
fraction of central galax;ies that are blue, one obtains an es- 
timate of the fraction of galaxies that turned red after they 
became satellite galax;ies (i.e. because of the environmen- 
tal quenching). An equivalent quantity can be defined using 
the measured fractions of pass ive and active satelli t e and 
central galaxies. In their study, Ivan den Bosch et al.l (|2008l ) 
find that satellite quenching affects roughly 40 per cent of 
all galaxies that are still blue at the time of accretion, inde- 
pendently of their stellar mass. 

Our results are shown in the top panel of Figure 1121 
and are compared to the same theoretical predictions that 
are given in Figure 1111 Top panels show the red transition 
fractions (Eq.[21), while bottom panels show the correspond- 
ing measurements for passive fractions. As the errors on the 
transition fractions are dominated by the errors on the frac- 
tion of red/passive satellites, we have shown here the same 
error bars used in Figure [TT] Both our red and passive transi- 
tion fractions increase as a function of halo mass. A similar 
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Figure 12. Lines show the same theoretical predictions given in Figure [TT] Data points with error bars show the red (top panels) and 
passive (bottom panels) transition fractions, i.e. an estimate of the fractions of satellites that have been quenched by the environment 
(see text for details). 



trend was recently found by IPeng et al.l (|201ll ) who used, 
however, local overdensity as a proxy for the environment. 
The red transition fraction also shows a clear dependency 
as a function of the galaxy stellar mass: for the two lowest 
mass bins considered, the red transition fractions increase 
from ~ 55 per cent in haloes of mass ~ 10^"^ M© to ~ 75 per 
cent in the most massive haloes considered in our sample 
(- 5 X 10" M0). For the most massive galaxies considered, 
no galaxy has been affected by environmental quenching in 
the lowest halo mass in our sample, while the transition frac- 
tion increases to ~ 60 per cent in the most massive haloes 
considered. The passive transition fractions are offset low 
with respect to the corresponding red fractions except for 
the highest stellar mass bin, and the dependency on stellar 
mass is weaker. On average, the passive transition fractions 
vary between ~ 30 per cent in the lowest mass haloes in 
our sample to ~ 65 per cent in the most massive haloes. For 
haloes more massive than 8 x 10^'' Mq , the estimated passive 
transition fractions follow quite nicely the trend expected for 
galaxies that spent more than ~ 7 Gyr in haloes more mas- 
sive than 10^^ Mq . For lower mass haloes, the figure favours 
a value of Thaio intermediate between 5 and 7 Gyr. 



5 DISCUSSION 

In the previous sections, we have characterized the envi- 
ronmental history of group and cluster galaxies. In partic- 
ular, we have related the observed fractions of passive/red 
galaxies to the fraction of galaxies that have spent more 
than a given time (Thaio) in a halo with mass larger than a 
given threshold (Mhaio). In this section, we discuss the main 
caveats of our analysis, and the main implications of our 
findings. 



5.1 Caveats 

As mentioned in Section^ the galaxy formation model used 
in this study over-predicts the num ber of low to inter mediate 
stellar mass galaxies. In particular. IGuo et al.l(|201ll ) showed 
that the model used here over-predicts the number of galax- 
ies in the stellar mass bin log(Mstar/Af0) = [9 — 10] by a 
factor ~ 2. For larger stellar masses, the model reproduces 
quite nicely the observed galaxy stellar mass function, as 
well as the halo occupation distribution measured for galaxy 
clusters both in the local Universe and at higher redshift 
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l|Poggianti et al.ll2010l . but see also iLiu et al.llioiol 'FI. If the 

excess measured in the stellar mass bin log(Mstar/M0) = 
[9 — 10] is due to galaxies that should have been stripped 
and/or merged with the central galaxy of their own parent 
halo, then there is likely an excess of galaxies with long sur- 
vival times as satellites. Since we base our analysis on the 
fraction of galaxies that spent more than a given time in 
haloes more massive than some threshold, our predictions 
for this mass range should not be significantly affected. 

We stress that our results are valid 'on average'. As 
discussed above, there is a relatively large halo-to-halo vari- 
ance that contributes, at least in part, to the scatter in the 
observed measurements. In addition, one should note that 
environmental effects very likely happen in a 'probabilistic 
way': not all satellites will continue forming stars for some 
given time and then simultaneously shut off. Otherwise, we 
would not observe a bimodality in the distribution of specific 
star formation rates. In reality, there will be satellite galax- 
ies whose star formation rate will be suppressed on relatively 
short time-scales and other satellites that will maintain some 
level of star formation activity for longer times. This is prob- 
ably in part related to the infalling orbital distribution of 
the satellite galaxies, as the orbital p arameters wil l influ - 
ence their dynamical friction timescale. iMcGee et al.l (|2009| ) 
argued that this should not affect significantly the results 
obtained if quenching requires relatively long time-scales, 
as those suggested by our study. The argument is based on 
the relatively narrow distribution of times that infalling dark 
matter subst ructures take t o reach t heir pericentr es from the 
virial radius (Ben sonll2005l . see also I Wetzel' 2011). It should 
be noted, however, that the evolution of galaxies infalling 
onto larger haloes depends on other factors than their or- 
bits, for example on the initial incli nation of the discs wit h 
respect to their orbital planes (e.g. IVillalobos et al.ll20l3 ). 
and variation in the inter-galactic medium properties of the 
host. Detailed numerical simulations are therefore needed to 
quantify the probability that a galaxy that has been orbit- 
ing within a halo of given mass has its star formation rate 
suppressed below some given threshold. As discussed above, 
our results suggest that the probability to be quenched ap- 
proaches unity for galaxies that have spent 5 — 7 Gyr in 
haloes more massive than 10^^ M©. These results might be 
useful in guiding the construction of halo occupation dis- 
tribution models based on abundance matching for passive 
and active galaxies. 

In Section U] we have m easured the red and passiv e 
transition fractions following Ivan den Bosch et al.l (|2008h . 
and assumed that the present-day population of central 
galaxies can be considered representative of the progen- 
itors of present da y satellite galaxies. As discussed in 
Ivan den Bosch et al.l (|2008l ). this assumption is subject to 
a number of caveats. In particular: (i) in reality, one should 
consider the population of central galaxies at the average 
time of accretion of the satellite galaxy population; and (ii) 
the stellar mass of satellite galaxies might increase if they 
continue forming stars after accretion, and might decrease 



^ Note that IPoggianti et al used a magnitude limited 

sample both in their observational and model catalogues. We have 
verified that their magnitude limit corresponds approximately to 
a stellar mass limit of ~ IQI^Mq. 



because of tidal stripping. Ivan den Bosch et al.l l|2008l ) ar- 
gue that the cumulative effect should be small. We stress, 
however, that satellite galaxies in our model have been ac- 
creted over a wide range of cosmic epochs (see top left panel 
of Figure O, and that the observed passive fractions require 
relatively long time-scales for the transition from blue/active 
to red/passive to occur. Therefore, the systematic errors on 
the transition fractions estimated in Section |3] could be rel- 
atively large. 

5.2 What is the role of pre-processing? 

In hierarchical cosmogonies, galaxy clusters are assembled 
through the accretion of lower mass haloes. Cluster galaxies 
might have resided for some time in 'groups', making prepro- 
cessing in thes e systems a potentially impo rtant phase for 
their evolution (IZabludoff fc Mulchaevll 19981 ). This has been 
discussed recently in two different studies that have used re- 
sults from A'^-body simula tions to analyse the accretion his- 
tories of cluster galaxies (|Berrier et al.l 1200 9: McGe e et al.l 
[iooi). The conclusions obtained by these studies are quite 
different: Berrier et al. estimate that ~ 10 per cent of the 
galaxies residing today in haloes of mass ~ 10" Mq were 
accreted onto the cluster potential when residing in haloes 
more massive than 10^'' M0 (see their figure 2), and argue 
that pre-processing in a group environment is "of secondary 
importance for setting cluster galaxy properties for most 
clusters". McGee et al. estimate, for clusters of similar mass, 
that about ~ 28 per cent of the cluster galaxies have been 
accreted from haloes more massive than 10^'^ (see their 
figure 2) . They also show that this fraction increases to ~ 45 
per cent for the most massive clusters in their sample, ar- 
guing that pre-processing in group environments "may be 
significant" . 

The reason for this discrepancy lies mainly in the adop- 
tion of two different even t s corr esponding to accretion onto 
the cluster. iBerrier et al.l (|2009l ) use A-body simulations of 
galaxy clusters (with resolution similar to that of the Mil- 
lennium Simulation used in this study), and assume that 
all bound substructures whose mass at the time of accre- 
tion is larger than some specific threshold, host a galaxy at 
their centre. By construction, the time used in Berrier et 
al. corresponds to our tsat, i-e. the time when a galaxy first 
becomes satellite of a larger halo. iMcGee et al.l (2009 ) use 
an approach that is equivalent to ours but base their analy- 
sis of a different semi-analytic model. Similarly to what we 
do, they trace the most massive progenitor of each galaxy 
back in time and record the halo mass of this progenitor at 
the time-step just before "the galaxy becomes a member of 
the final c luster", which shoul d correspond to our thaio. We 
note that iMcGee et al. I l|2009t ) did not distinguish between 
galaxies that are accreted onto the final cluster as centrals, 
and those that are accreted as satellites. Both studies as- 
sume that a cluster galaxy was 'pre-processed' if it spent 
some (significant) fraction of its lifetime in a 'group' halo 
with mass ~ 10^^ M© before 'being accreted onto the clus- 
ter', independently of whether the galaxy is a central or a 
satellite. 

As discussed above, the adoption of different charac- 
teristic times can lead to significantly different results, that 
also depend on the galaxy stellar mass. Considering only 
the haloes with mass ~ 10^* Mq in our sample, we find 
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that only ^ 1 per cent of the galaxies with stellar mass 
larger than 10^ M0 became satellites when residing in haloes 
more massive than ^ 10^^ Mq . The fraction increases to 
~ 4 per cent when considering galaxies more massive than 
~ 10^° M0 , and becomes ~ 37 per cent when considering 
the most massive galaxies in the sample (those with stellar 
mass larger than ~ 10^^ M©). Indeed, as we have discussed 
in Section [S] a larger fraction of these massive galaxies are 
accreted directly onto the final cluster as centrals. When 
considering the time when the galaxy is accreted onto the 
main progenitor of the cluster, we find that for the same 
haloes ~ 28, ~ 29, and ~ 44 per cent of galaxies with stel- 
lar mass larger than ~ 10^, ~ lO^'' and ~ 10^^ M0 were 
accreted from haloes more massive than ~ 10^^ Mq. Our re- 
sul ts are therefor e in very good agreement with those found 
by iMcGee et al] 12 009). An accurate comparison with re- 
sults from Berrier et al.l (|2009l ) is complicated by the use 
of a different approach that does not trace directly the 
evolution of the galaxy stellar mass. Their primary clus- 
ter sample contains about 16 galaxies per cluster, which in 
our case would correspond to a galaxy stellar mass limit of 
~ 2.6 X 10^" Mq. Considering the fractions estimated above, 
ou r results are therefor e not inconsistent with those obtained 
bv lBerrier eTa^\ (|2009l ). 

If we define as pre-processed all galaxies that have spent 
time as satellites of a lower mass system before becoming 
part of a clustei[f|, our results show that the fraction of galax- 
ies that can be pre-processed in a group-size halo of mass 
~ 10^^ Mq is significant. This fraction is largest for lowest 
mass galaxies. As discussed above, while ~ 44 per cent of the 
most massive galaxies have been accreted from haloes more 
massive than ~ 10^"^ Mq , a large fraction of those have been 
accreted as centrals. On the other hand, basically all the 
~ 28 per cent low and intermediate mass galaxies that have 
been accreted from haloes more massive than ~ 10^^ Mq are 
satellites at the time of accretion. 

We stress that these results are natural consequences 
of hierarchical structure formation, and that they do not 
depend significantly on the particular semi-analytic model 
that we have used in our study. In order to understand 
how this preprocessing can affect the observed morphology- 
density relation, morphological mix, etc. it is therefore cru- 
cial to study the relevance and effects of various physical 
processes, at the scales typical of galaxy groups (specifi- 
cally, halo masses ~ 10^'^Mq). Published numerical work 
on the role of environment on galaxy evolution, however, 
has so far focused largely on massive galaxy clusters, with 
very few numerical studies devoted to understand the effect 
of the group environment. Further numerical work in th is 
direction is clearly needed (see e.g. IVillalobos et al.ll20l3 ). 



^ We stress that this definit ion differs from that adopted in 
some earlier studies, including iBerrier et al. I ||2009|' I iMcGee' et al.l 
who did not distinguish between physical processes acting 
only on satellite galaxies and processes that can affect all galaxies 
in a group. 



5.3 Environment and time-scale of galaxy 
transformations 

It has long been known that the local galaxy population con- 
sists roughly of two different types of galaxies: red galaxies 
with low levels of ongoing star formation and blue galax- 
ies with active star formation. This bimod ality is known to 
extend at least up t o z ~ 1 or higher (|Bell et al.l |2004| : 
Whitakcr ct al. 2011), and to depend on the environment 
(e.g. Baldry et al. ,2006 ) ■ The physical origin of the ob- 
served bimodal distribution remains a question to be an- 
swered. In particular, it is unclear what is the characteristic 
time-scale of the galaxy transformation from blue/active to 
red/passive, and if this transformation is associated with a 
particular environment. 

In our study , we have used the group catalogue by 
lYang et al. ( 20071) complem ented with a cluster catalogue 
(*von der Linden et al.l [20071 ) . and have tried to constrain 
these processes using two observational constraints: how the 
fraction of passive galaxies varies as a function of halo mass, 
and how the same quantity varies as a function of cluster- 
centric distance. In order to estimate the efficiency of satel- 
lite quenching, we have used the same statistics discussed by 
Ivan den Bosch et al.l l|2008l ). We have shown that the frac- 
tion of satellite galaxies that have been quenched by the 
environment increases as a function of halo mass, from ~ 30 
per cent in the lowest mass haloes in our sample to about 
~ 65 per cent in the most massive haloes considered. The 
dependency on stellar mass is not large, and it tends to de- 
crease for increasing halo mass. 

Our calculations show that these observational trends 
are naturally explained by the fact that a larger fraction of 
satellite galaxies have spent long times in relatively massive 
haloes. Therefore, a more efficient quenching in more mas- 
sive haloes is not required. These considerations give little 
support to cluster-specific processes like ram-pressure as the 
main drivers for the observed trends. In addition, the com- 
parison with our theoretical predictions show that satellite 
galaxies become passive after they have spent a time Thaio ~ 
5 — 7 Gyr in haloes more massive than Afhaio ^ 10^'' Mq . 

Several previous studies have argued for rela- 
tively long timescales for s atellite quenching. E.g. 
iBalogh. Navarro fc Morr is (2000) study the origin of cluster- 
centric gradients in star formation rates and colours of rich 
clusters assuming they are built from the accretion of field 
galaxies, and argue for a gradual decline (over a timescale 
of a few Gyr) o f the star formation rate of cluster galaxies 
after accretion. IWang et al.l (|2007h use a physically based 
halo occupation distribution model and the observed distri- 
bution of 4000 A break strengths to constrain the star for- 
mation histories of central and satellite galaxies. They find 
that satellite galaxies have declining star formation rates, 
with average e-folding times Ts ~ 2.5 Gyr, and no signif- 
icant dependency on stellar mass. iKang fc van den BoschI 
l|2008l ) use a semi-analytic model of galaxy formation and 
the observed increase in the red fraction of satellite galax- 
ies as a function of stellar mass. They argue that model 
results can be brought into good agreement with the ob- 
servational data by decreasing the stripping efficiency of 
the hot gas reservoirs associated with galaxies being ac- 
creted onto a larger halo. In particular, they assume a con- 
stant stripping rate over a time-scale of ~ 3 Gyr. Recently, 
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IWetzel. Tinker fc Conrwl HoTJ) have used a group cata- 
logue based on SDSS DR7 and studied the distribution of 
SSFRs for satelhte galaxies and its dependence on stellar 
mass, halo mass and halo-centric distance. The 'persistent 
bimodality' they find for satellite galaxies indicate that star 
formation in active satellites continues to evolve (as in ac- 
tive centrals) for several Gyr. Several other studies have ar- 
gued for relatively long timescales for the suppression of the 
star formation rates in satellite galaxies jFinn et alJlioO^: 
Weinmann et al.l l2009l: iDe L ucia et alj |2009|: I Simard et alj 



20091 : 1 von der Linden et al.| l2010: Mc Gee et al.ll201ll ') 



As explained in Section [2] a gradual stripping of the 
hot gas reservoir has been advocated and assumed in 
many recent studies based on semi-analytic models (see e.g. 
iFont et al.l[200i : I Weinmann et al.ll201ol : Iguo et al.ll201lh . It 
is, however, unclear how such a gentle mode of strangulation 
could support cooling and star formation in satellite galaxies 
for several Gyr. More detailed numerical studies are desir- 
able to understand if the main problem with satellite galax- 
ies in semi-analytic models is related to a poor treatment of 
environmental processes, or if it signals a more fundamental 
problem with the treatment of star-formation and feedback 
for these galaxies. 



6 CONCLUSIONS 

We exploit publicly available catalogues from semi-analytic 
models to study the environmental history of group and 
cluster galaxies. In particular, we focused on haloes with 
M200 > 10^^ M© , carried out our analysis within a theoret- 
ical framework, using the halo mass (M200) as a proxy for 
the environment, and considered explicitly the dependency 
on galaxy stellar mass. Our analysis highlights a number 
of natural consequences of structure formation that are im- 
portant to consider when interpreting observational data. In 
particular, we show that: 

• On average, the surviving massive satellites within 
galaxy groups/clusters were accreted later than their less 
massive counterparts and they come from more massive 
haloes. The most massive galaxies tend to be accreted 
onto the main progenitor of their final group/cluster when 
they are central galaxies of their own haloes. Less massive 
group/cluster members become satellites over a wide range 
of redshifts, and about half of them are accreted onto the 
final group/cluster when they are already satellite galaxies. 

• The mixing of galaxy populations is incomplete during 
cluster assembly, which establishes a correlation between the 
time a galaxy is accreted onto a more massive halo (i.e. it 
becomes a satellite galaxy) and its distance from the final 
cluster centre. The radial trend is weakest for the most mas- 
sive galaxies because of efficient dynamical friction and the 
late formation times of massive haloes. 

These trends can be considered as the result of a his- 
tory bias, that represents an integral part of the hierarchical 
framework. According to the current paradigm for structure 
formation, dark matter collapses into haloes in a bottom-up 
fashion. Small systems form first and subsequently merge to 
form progressively larger systems. As structure grows, galax- 
ies join more and more massive systems, therefore experienc- 
ing a variety of environments during their lifetime. In this 



framework, where a galaxy resides today (or at any cosmic 
epoch) depends heavily on what its past history was. Our 
analysis demonstrates also that binning galaxies according 
to their stellar mass does not suffice to disentangle the role 
of nature and nurture because galaxies of different stellar 
mass have different environmental histories. So, for exam- 
ple, two galaxies of identical mass at some cosmic epoch can 
end up having different stellar masses if one of them falls 
onto a cluster and the other remains in a region of average 
density. 

The traditional nature versus nurture debate, as well as 
the controversy regarding the primacy of stellar mass versus 
environment, represent then subtle issues regarding 'corre- 
lations' and 'attribution'. It is, in principle, possible to sep- 
arate nature vs nurture if they are correlated but physically 
uncoupled. However, our results demonstrate that the two 
are strongly and physically connected so that any attempt 
to separate them is ill posed and misguided. We stress that 
this does not mean that studying galaxy properties at fixed 
stellar mass but in different environments does not provide 
important information, as demonstrated by a number of re- 
cent studies, many of which are referred to in this work. 
Results, however, should be interpreted with care. 

To understand and quantify how hierarchical structure 
formation affects the observed environmental trends, we 
considered two observational constraints: the dependence 
of the passive galaxy fraction on halo mass, and cluster- 
centric distance. We based our observational measurements 

and the clus- 



on the group cat alogue bv | Yang et al.| ( 200 
ter catalogue bv lvon der Linden et al.l (200 
to-date estimates of the star formation 



We used up- 



rate and stellar 
masses, and complemented the cluster data with observa- 
tional data availab le for the Com a and Virgo cluster pre- 
sented by Weinma nn et al.l |20l3). 

We find that significant fractions of group / cluster galax- 
ies have been accreted onto their final halo as members of 
groups with mass ^ 10^'' Mq . For the particular model used 
in our study, and considering all haloes in our sample, we 
find that this fraction is ~ 27 per cent for galaxies with stel- 
lar mass ~ 10^ and ~ 10^° Mq , and increases to ~ 44 per 
cent for galaxies with mass ~ 10^^ M©. We also find that 
~ 48, ~ 43, and ~ 23 per cent of the galaxies in the same 
stellar mass bins are accreted onto the main progenitor of 
the final halo as satellite galaxies. Large fractions of group 
and cluster galaxies have therefore been 'pre-processed' as 
satellites of groups with mass ~ 10^'^ M©. Since a large frac- 
tion of the most massive galaxies we see in clusters today 
have been accreted as centrals, they have been pre-processed 
the least. 

Comparisons with observational data suggest that 
galaxies become passive only after they have spent Thaio ~ 
5 — 7 Gyr in haloes more massive than Mhaio = 10^^ M©. The 
observational trends are naturally explained by the growth 
of the large scale structure with no need for a more effi- 
cient quenching in the most massive haloes. It is unclear how 
satellite galaxies can sustain significant levels of star forma- 
tion for such long time-scales and if this can be achieved, 
within the framework of current semi-analytic models of 
galaxy formation, by simply relaxing the assumption of in- 
stantaneous stripping of the hot gas reservoir associated 
with galaxies when they are accreted. On the numerical side, 
very little work has been carried out to understand the im- 
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portance and effects of different physical processes at the 
velocity dispersions typical of the group environment. Fur- 
ther studies in this directions are needed. 

In this study, we have limited the analysis to haloes 
selected at redshift zero. A more comprehensive analysis ex- 
tending to higher redshifts is clearly desirable and will likely 
provide important constraints on the physical processes (and 
related time-scales) responsible for establishing the observed 
environmental trends. Finally, we stress that our study is 
based on theoretical definitions of environment. In future 
work, we plan to extend the analysis outside the virial radius 
of dark matter haloes, and to use environment definitions 
that are closer to those commonly adopted in observational 
studies (e.g. local density). We believe that such an anal- 
ysis will provide important guidance in the interpretation 
of observational results in a cosmological context, and in 
the comparison of studies at different cosmic epochs and/or 
from different surveys. 
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